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Abstract 

We consider pulses of finite duration for coherent control in the presence of classical noise. We derive the 
Q_| corrections to ideal, instantaneous pulses for the case of general decoherence (spin-spin relaxation and spin- 

lattice relaxation) up to and including the third order in the duration Tp of the pulses. For pure dephasing 
(spin-spin relaxation only), we design tt and 7r/2 pulses with amplitude and/or frequency modulation which 
resemble the ideal ones up to and including the second order in Tp. For completely general decoherence 
including spin-lattice relaxation the corrections arc computed up to and including the second order in Tp as 
well. Frequency modulated pulses arc determined which resemble the ideal ones. They are used to design 
a low-amplitude replacement for XY8 cycles. In comparison with pulses designed to compensate quantum 
noise less conditions have to be fulfilled. Consequently, we find that the classical pulses can be weaker and 
simpler than the corresponding pulses in the quantum case. 



'Corresponding authors 
Email addresses: christopher.stihlStu-dortmuiid.de 
(Christopher Stihl), benedikt . f ausewehOtu-dortmund. de 
(Benedikt Fausewoh), s .pasiniOf z-juelich. de (Stefano 
Pasini,), goetz.uhrigatu-dortmund.de (Gotz S. Uhrig,) 



Preprint submitted to Journal of Magnetic Resonance 



October 17, 2012 



Modulated pulses compensating 
classical noise 



Christopher Stihl^, Bcnedikt Fauscwch'^, Stcfano 
Pasini,'^'*, Gotz S. Uhrig,'^^* 

Lehrstuhl fiir Theoretische Physik I, TU Dortmund, 
Otto-Hahn Strafie 4, 44'^'^^ Dortmund, Germany 
Forschungszentrum Jiilich, 52425 Jiilich, Germany 



1. Introduction 

In open quantum systems the interaction with an 
environment is one of the main reasons for the decay 
of the coherence and for the loss of information, both 
obstacles to an efficient processing of quantum infor- 
mation or to high precision NMR measurements. But 
by means of a proper time-dependent modulation of 
the system dynamics, it is possible to average out the 
effects of the environment. This leads to the prolon- 
gation of the coherence time and to the suppression 
of the dccohcrence. 

In NMR, a standard way of time-dependent mod- 
ulation is the application of very short pulses which 
flip spins. This approach dates back to the spin echo 
of Erwin Hahn in the fifties [l[ and then was quickly 
extended to periodically applied pulses by Carr and 
Purcell and Meiboom and Gill Q. Since then 
a large variety of periodic pulse sequences has been 
proposed 0, [Ij among which we highlight the XY4 
cycle 0, 0] and its derivatives as a means to preserve 
all three spin directions. 

Recently, the approach of stroboscopic pulsing of 
two-level systems to suppress their decoherence has 
been revived under the name of dynamic decoupling 
(DD) for open quantum systems in the context of 
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quantum information processing (QIP) d, H, [lo| . 
Again, the idea is to effectively decouple the two-level 
system from the environment which induces its de- 
coherence. Interestingly, quantum information pro- 
cessing requires a slightly differing performance than 
NMR. During the information processing it is desired 
that the decoherence is much better suppressed to the 
level of 10~^ to 10~^. But deviations at longer times 
do not matter. In contrast, in NMR a loss of signal of 
the order of 50% still allows one to track the signal, 
but it is desired to be able to do this for as long as 
possible. 

Hence, in QIP novel sequences of pulses have been 
recently developed which do not rely on periodically 
repeated cycles, but vary the time intervals between 
the pulses. A first suggestion relies on recursively 
concatenated sequences (CDD). They require an ex- 
ponentially growing number of pulses if longer 
times have to be reached. For the suppression of 
dephasing (spin-spin relaxation) the Uhrig dynamic 
decoupling (UDD) is much more efficient because it 
requires oiily a linearly growing number of pulses 
fl4 13, 3 [ill- practice, this pays for power spec- 

Subsequently, UDD 
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has been extended to also suppress general decoher- 
ence including spin-lattice relaxation due to spin flips. 
This sequence goes under the name of quadratic dy- 



namic decoupling (QDD) [18|, llE l20|, l2ll and it is 
further generalized to nested UDD to deal with more 
than one two-level system, i.e., with quantum regis- 
ters [22| . Experimental verification of these theoreti- 
cal proposals are also available for the UDD [23. 24 1. 

All the above approaches rely theoretically on 
ideal, instantaneous pulses which can not be realized 
experimentally. There are various routes to overcome 
this problem. One is to abandon the use of pulses al- 
together by resorting to continuous control fields, see 



for instance Ref. [25|. We will not follow this route 



here but stick to sequences of pulses. We point out, 
however, that there is a crossover between both ap- 
proaches if the pulses become so long that one touches 
the next one or if a pulse is designed such that it acts 
itself like a cycle of pulses, see below. 

For pulses, the most direct solution is to realize 
very short pulses so that their duration Tp is so short 
that the coupling to the bath does not play a role 
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during the pulse 



;7p(Tp,0) = P,^+O(rp) 



(1) 



where Up (Tp,0) is the unitary time evolution under 
the combined action of the coupling between system 
(spin, qubit, or two-level system) and bath as well 
as the coherent control while stands for the time 
evolution under the pulse alone as if the system were 
isolated. 

A more sophisticated solution consists in shaping 
the pulse such that the above approximation holds in 
the form 



[/p(Tp,0) = P.,,+O(rp™ 



(2) 



We will call a pulse which fulfills ^ an m-th order 
pulse. Thus, a standard unshaped pulse is a 0-th 
order pulse. In this work, we will present first and 
second order pulses. 

We emphasize that pulses with the property ^ 
cannot be used as simple replacements for an ideal, 
instantaneous pulse because they still last for the 
time Tp. A proper replacement would have to fulfill 

f/p (Tp, 0) = U (Tp, Tid) Pr,U (Tid, 0) + 0(t™+1) (3) 

where U (t2, ti) is the unitary time evolution without 
any external control and Tid is the instant at which 
the 5-like control pulse fV^ is inserted. The time 
evolutions U before and after an ideal, instantaneous 
pulse Pj-p simulate that right before and after an ideal 
pulse the system is coupled to its decohering environ- 
ment. It could be shown, however, that for tt pulses 
the ansatz ([3]) can at most be fulfilled for m < 1 so 
that the route of findin g pr oper replacements cannot 
be pursued further 26l. 27|. 

Hence, the ansatz ^ had to be developed further 
because it can be realized for arbitrary m in principle 
[28| and concrete proposals have recently been made 
for second order pulses (to = 2) for general quantum 
baths [2^,[33l- I* could be shown that pulses with this 
property can be used in sequences of pulses if they 
are not used in the same manner as instantaneous 
pulses would be used, but with adapted sequences 
[3ll [32 1 . These observations clearly show that the 



Of course, shaping of pulses has also been a 
long-standing issue in NMR. The main aim was to 
generate robust pulses which are only weakly sus- 
ceptible to pulse imperfections. We can mention 
only partly the abundant literature on this issue 

Mm, [Mjii, [n m, m, HO, im, m m, for a book 



see Ref. [44 1 



In contrast to earlier work, e.g., the one by Skinner 



shaping of pulses is an important ingredient for high- 
fidelity coherent control. 



and collaborators [37|, |38[ , we focus in this article on 
analytical derivations as far as possible. In contrast 
to our previous work on pulses for systems coupled to 
quantum baths 2^, we will discuss classical baths. 
The fundamental reason is that in most experiments 
the dominating fluctuations destroying coherence are 
of classical nature, see, e.g., Ref. [45|. Generically the 
decohering fluctuations are induced by a large num- 
ber of microscopic and macroscopic degrees of free- 
dom. These degrees of freedom are at least partly 
at rather high temperatures relative to their generic 
energy scales. For instance, the nuclear spins can 
mostly be considered to be in a disorderd state corre- 
sponding to infinite temperature. Thus the resulting 
fluctuations are thermal fluctuations and their quan- 
tum character plays only a smaller role. 

One may think that the strongly disordered states 
of the decohering baths pose a problem to the preser- 
vation of coherence in the systems under study. But 
the opposite is true: The lack of quantumness of the 
fluctuations allows us to consider classical fluctua- 
tions only. This implies less restrictive conditions on 
TO-th order pulses of the type ([2|). Hence, the neces- 
sary pulse shapes are simpler and most importantly 
the required amplitudes are lower than in the full 
quantum case. The main goal of the present article is 
to show this explicitly. The concomitant message to 
experiment is that for given power supply the classi- 
cal pulses will be easier to realize. They can be made 
shorter than their quantum analogues. 

In this article, we derive amplitude and frequency 
modulated pulses for pure dephasing due to classi- 
cal noise, i.e., for spin-spin relaxation without spin 
flips, up to and including second order (m = 2). Par- 
ticular attention will be paid to the minimization of 
the amplitude because this is important for practical 
use. For general decoherence due to classical noise, 
i.e., including spin flips due to spin-lattice relaxation. 
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we derive frequency modulated pulses up to and in- 
cluding second order (m = 2). We stress that such 
pulses can be used to build continuous versions of the 
well-known XY8 cycle [1,01. 

The paper is organized as follows: In section [2] the 
equations for the decoupling of general decoherence 
are derived. In sections |3] and H] these equations are 
shown for the first and second order. They are solved 
numerically for the case of pure dephasing, both for 
amplitude and for frequency modulated pulses. In 
section [5] we study decoupling pulses that combine 
the characteristics of amplitude (AM) and frequency 
modulation (FM) in order to illustrate that transient 
amplitudes, when the pulse is switched on or off, can 
be easily dealt with. In section [S] the case of gen- 
eral decoherence is solved for by frequency modulated 
pulses. We draw our conclusions in section [71 

2. General Equations 

2.1. Ansatz 
We consider the following Hamiltonian 

Htot (t) =H{t)+ Ho {t) (4) 

which consists of two parts: A system Hamiltonian 

H{t) = rJ{t)-<j (5) 

and a control Hamiltonian 

Hoit)^v{t)-a. (6) 

The system Hamiltonian describes the physics of 
a single spin coupled to a classical bath through the 
time-dependent random function f]{t). We describe 
decoherence by averaging over f](t) . One may assume 
f]{t) to represent Gaussian fluctuations, but this is 
not essential for the present work. The vector a with 
components given by the Pauli matrices ax, cFy and 
az stands for the spin operator. The spin is subject 
to general dephasing including spin-lattice relaxation 
if all three vector components are present. 

In the control Hamiltonian Holt), v(t) is the time- 
dependent vector- valued amplitude of the pulse op- 
erating on the spin. At each instant in time, the spin 



undergoes a rotation about the time-dependent axis 

In all what follows we assume that the control \v\ 
is large enough tt « vTp so that Tp can be chosen 
small enough and tt ^ ||i?o||Tp holds. Then it is well 
justified to expand in the decoherence dynamics. For- 
mally, we perform an expansion in Tp while |u|Tp is 
kept at order unity. 

Following the technique developed in Refs. [2^, 
[30j we split the complete time-evolution operator 

Up (Tp,0) = 7~|cxp[— i / Htot (^dt] I of the system 

during a pulse of duration Tp into two terms: (i) The 
evolution of the spin under the effect of only the con- 
trol field Pr - as if it were decoupled from the bath 
- and (ii) the correction Uc (rp, 0) 

f/p(Tp,0) :=F,^[/e(Tp,0). (7) 

The notation T stands for the quantum mechanical 
time-ordering operator. The operator J7c(Tp,0) in- 
corporates the deviations from an ideal rotation re- 
sulting from the interaction between the spin and the 
bath. 

The rotation operator can be expressed as 

A :-cxp{-i<T.a(t)^}, (8) 

where the unit vector a stands for the axis of to- 
tal rotation while ip{t) is the corresponding angle by 
which the spin is rotated from the instant t — till 
the instant t ^ t. 

Both operators Pt and Up {t, 0) satisfy a 
Schrodinger equation at every instant t, for de- 
tails see Ref. [2^ [s^. For the former operator 
idtPt = Ho{t)Pt translates to 

2v(t) = 4>'{t)a{t) -f a'{t) sin'0(i) 

- [l- cos at)) [a' {t)xa{t)]. (9) 

For the latter, one obtains the formal solution for 

i e [o,Tp] 

C/e(rp,0)=r{e-'/o^''^W'i*}, (10) 

where 

H{t):=P,-^H{t)Pt. (11) 
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2.2. Integral Equations 

One of the standard procedures to deal with time- 
dependent Hamilton operators as in (|10p is the Mag- 
nus expansion 46|, |47[. It consists in replacing the 



time dependent Hamiltonian by a sum of average 
time-independent Hamiltonian terms 



U (0, Tp) = exp 




(12) 



where = 0(Tp ^). The explicit form of the first 
three cumulants of the expansion is given by 



(13a) 

H{tx),H{t2) (13b) 
dt. 









Tp#i = 


r 

Jo 


dtH{t) 




i 




Tp^2 = 




fdti 
Jo 








Tp^3 = 


"6 


fdti 
Jo 



{[[H{t,),H{t2)] ,Hit3) 
+ [H{t,),[H{t2),H{t3)W} 



(13c) 



The rotated spin operator PtdPt can also be un- 
derstood as the vector a rotated by the angle tp 
about the axis a{t). By applying this rotation ma- 
trix not to a, but its inverse to fj we can reex- 
press the Hamiltonian as (t) = Hfj (t) • tr, where 
riff (t) := [Dg {—ip{t)) ■ r]\. The corresponding 3 x 3 ro- 
tation matrix Dg [tji) is given in [Appendix B[ The ro- 
tated vector 77 has the following explicit form 2^, S^] 



it) = cos (1/') rf — sin (?/») [a x ^ 
+ [1 — cos (V')] {ff ■ a)a 



(14) 



where we did not denote the time dependence of -0, 1^^ 
and a for the sake of clarity. 

Now we are able to derive an expression for the 
terms of the Magnus expansion 



TpHi = y dt nff{t) ■ (T, 




(15a) 



TpH2 = / dti / dt2 a{n,j{ti) x n^(i2))- (15b) 



ti 



dti / dt2 / dtg 



Clearly, if 



[2nff{t2) {riff {ti)fiff (t^)) 
~ nffits) {fiif (ti) n,j {t2)) 
- n,j{ti) {fiij (tz) n,j {t2))] 



Hi=H2=Hs^ 0, 



(15c) 

(16) 
3. 



holds, the corresponding pulse fulfills ([2]) for m 
If only Hi and H2 vanish then ([2]) holds for m = 2. 

But we do not know the time-dependent ri(t). As 
usual in statistical physics we average over all possible 
realizations {ri(t)} for which we need the correlations 
of {r]{t)}. But it is in general not sufficient to average 
Uc over {ri{t)}. One has to average the changes in 
the density matrix p of the quantum spin. Thus we 
consider U^pUp and average this quantity over {ri{t)} 



UlpUp {uIpUv)v 



(17) 



as indicated by (■),,. For S = 1/2, i.e., the two-level 
quantum system under study, p can be expanded in 
the Pauli matrices and the identity. The latter is 
transformed trivially so that it is sufRcient to consider 

{UlaUp)r, = {Ula'U,)r, (18) 

where a' = P^^aP-r^. To order m we require that 

(C/ta'[/e), = a' + 0(r™+i). (19) 

The vector a parametrizes all deviations of the den- 
sity matrix from the identity and so does a' . Hence 
we need not distinguish a and a' . 

In an expansion up to second order in Tp we verified 
that our objective ([T^ is indeed fulfilled if {Uc),-i is 
the identity up to terms of the order ©(t™^^). Thus 
in the sequel, we will consider the average of the uni- 
tary {Uc)n- In third order, however, one would have 
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to consider ([TO)) rather than the average of Uc alone. 
Next, we will address the solutions for the conditions 
(Hi)^ = 0. To do so we will specify these conditions 
further. 

3. Amplitude-Modulated Pulses for Pure De- 
phasing 

Here we consider pure dephasing which means that 
we exclude spin flips. Thus we restrict the Hamilto- 
nian of the system to 



H{t) = rjit)a. 



(20) 



In this case rotations around the y axis in spin space 
(a = (0, 1, 0)) are sufficient to decouple the spin from 
the bath. Thus the control Hamiltonian may assume 
the following simple form 



Ho{t)=vit) ay 



(21) 



so that we only have to consider the amplitude mod- 
ulation of v(t). Despite its simplicity this approxima- 
tion is representative for a large class of the decoher- 
ing systems in which Ti is much larger than T2. This 
can be reached by large magnetic fields which split 
the two energy levels strongly. Then all terms differ- 
ent from are averaged out in the rotating frame 
approximation. 

For Hq as in (PT|) the rotation operator Pf = 
cos [-0(0/2] — i(Ty sm[ip{t)/2] transforms the system 
Hamiltonian (gOl) to H{t) in ([TTl) yielding 

H{t) = 77(0 [cos{i;it))a, - sin(0(O)(T,] . (22) 

From Eqs. ([Ti]) . we obtain 

^vit) = V{t) [" sin(V'(0), 0, cos(V'(i))] ■ (23) 

Due to the classical nature of the bath only a few 
of its average values are needed. In the leading order 
only the mean value 



(24) 



enters. Since we assume the noise to be uniform in 
time, i.e., it does not depend on time, the mean value 
is a constant. We find that the vanishing of TpHi = 



fJ-i,i'^x + requires the vanishing of fiij with 

i £ {1,2} with 

V sin(-i/'(t)), (25a) 

Jo 

Aii,2 :=?7 / "dt cos(V(0), (25b) 
Jo 

unless fj is zero accidentally. As expected the first- 
order corrections are of order Tp. 

It is straightforward to find solutions of ([?51) for 
piecewise constant amplitudes or continously vary- 
ing amplitudes. Since the first order equations (|25|) 
arc identical to the ones for the quantum case, we 
studied and solved them before [2^ . We refrain from 
presenting such solutions here again. We stress that 
solutions of the first order equations ([^5)) had been 
found before in the context of searching for pulses 
which are robust against frequency offsets 33,35, 3^ . 



For the second order in Tp we need some infor- 
mation from the autocorrelation function g{A.t) :~ 
{r}{t)ri{t + ^t))rj. For short delays Ai we can 
expand g{At) in the form 



g{At) = if + s^) + gi\At\ + 0{Aty^ 



(26) 



where s is the usual variance = ((77 — rff')^^ We 
do not exclude the appearance of the term propor- 
tional to I At I although it is at odds with analyticity. 
Often, a very fast microscopic process makes 77 (i) be- 
have as in an Ornstein-Uhlenbeck process which is 
characterized by a cusp in g[At) and a Lorentzian 
power spectrum after Fourier transform, see for in- 
stance Ref. 48|. But we will see below that for the 
purposes of the present article even gi does not mat- 
ter for second order pulses. It occurs only in third 
order pulses. 

Explicit calculation yields the second order correc- 
tion 

TpJ?2 ■^cTyiT'f + s^)- 



dti / dhsm[iih)-i;{t2)], (27) 
Jo 



which is quadratic in Tp. 

Comparing the integral equations for the first and 
the second order with those derived for a quantum 
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Figure 1: (Color online) Symmetric second order tt pulse with 
piecewise constant amplitude suppressing pure dephasing. The 
symmetric switching instants are ri = 0.07623078rp , T2 = 
0.26784319tj,, = 1 — T2 and = 1 — t\, the amplitude is 



±6.72572865r"^ 



Figure 2: (Color online) Symmetric second order ^ pulse with 
piecewise constant amplitude suppressing pure dephasing. The 
symmetric switching instants are t\ = 0.03312609Tp, T2 = 
0.25209296rp, T3 = 1 — T2 and T4 = 1 — ti, the amplitude is 
±6.32709469t"^ 



bath [i^ we find that they are identical in first or- 
der. In second order, they are similar in form but 
in the quantum case there are three equations while 
two of them vanish in the classical case because the 
system Hamiltonian ((20|) commutes with itself at dif- 
ferent times. This implies that it is less demanding to 
find a numerical solution for the equations \2b\ and 
(|27|) in the classical case than in the quantum case. 
Practically, pulses can be found with a lower ampli- 
tude than in the quantum case. In Figs. [1] and [2] 
examples of a tt and a 7r/2 second-order pulses with 
piecewise constant amplitudes arc shown and com- 
pared to the corresponding pulses for a quantum bath 
[29[ . Clearly, the maximal amplitudes of the classi- 
cal pulses are lower. A technical remark is in or- 
der: Although the classical pulses have to fulfill less 
conditions we designed them with the same number 
of switching instants as the quantum pulses. But 
the absolute value of the amplitudes of the classical 
pulses is always the same while it varies in the quan- 
tum case. So the amplitude values arc the additional 
variables needed in the quantum case. 



represented by 



v{t) 



27r 
cos I 1 



'47r 

[b — a) cos ( — t 



ocos 1 



Equations (|25l27p can also be solved for continuous 
pulses. For symmetric pulses the function v{t) can be 



(28) 



where 6 is either tt or 7r/2 and a and h are constants. 
The pulse fulfills the requirements w(0) = w(Tp) = 
and ^'(O) — w'(Tp) = 0, as shown in Figs. ([3]) and 
(U]). Also in this case the maximal amplitude is lower 
than in the quantum case. 

All solutions presented above were obtained nu- 
merically. The piecewise constant solutions for am- 
plitude modulation were found using "fsolve()" from 
the "scipy" library for Python, which essentially is 
a wrapper around the "hybrd" and "hybrdj" algo- 
rithms from MINPACK. No further numerical cal- 
culations were needed to obtain those solutions be- 
cause the integrations were done analytically. This 
was achieved by partitioning the integration domain 
into the intervals between the switching instants. In 
this way ^{t) becomes a linear function within each 
interval of the partition such that the contribution of 
each interval is analytically available. 

To obtain continuous solutions for amplitude mod- 



7 




Figure 3: (Color online) Symmetric tt pulse with continuosly 
modulated amplitude correcting up to second order suppress- 
ing pure dephasing. The function v(t) is parametrized by Eq. 
12511 with a = -1.92179255 and b = 2.86838351. 



Figure 4: (Color online) Symmetric ^ pulse with continuously 
modulated amplitude correcting up to second order suppress- 
ing pure dephasing. The function v{t) is parametrized by Eq. 
l[28)l with a = -5.41258549 and b = -3.48909926. 



ulation, the GNU Scientific Library (GSL) was used 
for performance reasons. To find the multidimen- 
sional roots of the necessary sets of equations, the 
"gsLmultiroot_fsolver_hybrids" algorithm was used. 
It is very similar to the "hybrd" algorithm in MIN- 
PACK except for the internal scaling. A proposed 
solution is accepted if the residue, i.e., the sum over 
all absolute values is less than 10~^°. The 

integrand values fi themselves are calculated using 
"gsLintegration_gaq" which relies on an adaptive nu- 
merical integration using the 61 point Gauss-Kronrod 
rules until the estimate of the absolute global error is 
less than 10"^'^. To compute multiple integrals, calls 
of "gsl jntegration_gaq" were nested appropriately. 

4. Frequency Modulated Pulses for Pure De- 
phasing 

It is important to analyze the case of frequency 
modulated pulses as well because there may be ex- 
perimental setups where frequency modulation (FM) 
can be implemented much easier or more accurately 
than amplitude modulation. There is only an ini- 
tial and a final jump in the amplitude, the remain- 
ing control evolves smoothly. The initial and final 
transients will be discussed in the following section. 



Another advantage of frequency modulated pulses is 
that rotations about different axes can be realized in 
a natural way. Thereby, all kinds of deviations from 
the initial spin state can be compensated including 
spin relaxation. This latter point will be exploited in 
more detail in section [6l 

We study pulses acting only in the xy-plane of the 
spin orientation with a fixed amplitude Vq := \v\. In 
the rotating frame taking the Larmor frequency wl 
into account, the current axis of rotation is given by 

/ Vocos{n{t)) \ 
vit) = I Vo shi{n{t)) 1 (29) 

where n{t) is a time dependent phase which is tuned 
externally. The control Hamiltonian Hq is realized by 
applying a field perpendicular to the cr^-axis which is 
rotating with the Larmor frequency in the laboratory 
frame 

Hq""^ = Vq [a^ cos(wLi - ^(t)) - <7y sin(tJLi - ^{t))] ■ 

(30) 

The derivative dtil.{t) of the time dependent phase 
represents the deviation of the current frequency from 
the Larmor frequency. In this sense Eq. (|30p describes 
a frequency modulated pulse. In the rotating frame 
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we have 



Ho{t) ^ v{t) ■ a. 



(31) 



In order to find a{t) and ip{t) appearing in the 
paramctrization ^ of the pulse one has to solve the 
differential equation Because a{t) is a unit vec- 
tor, it can be parametrized by two angles (j){t) and 
9{t) 



a{t) := 




sm{e{t))cos{(j){t)) 
sin(6'(t)) sm{(j){t)) 
cos(6'(<)) 



(32) 



Solving Eq. ([9]) for the time derivatives of ij{t), (f>{t), 
and d{t) we finally find 



df ip = 2Vb sin 9 [sin ft sin (f> + cos O cos (j>] (33a) 



cos ^ sin(17 — 0) — sin ^ cos cos(ri 



sin Y sin ( 



(33b) 



dt9 = 



cos ^ cos 9 cos(51 ~ (j)) + sin ^ sin(17 



smf 



(33c) 



The seeming singularities for 9 ~ mn, to e Z on 
the right hand sides of Eqs. (j33bp and p3cp have no 
physical reason. They only result from the choice 
of spherical coordinates where (j) is ill-defined for 
9 ~ mn^m G Z. In contrast, the global axis of rota- 
tion a is ill-defined if i/j is a multiple of 27r because 
then the unitary of the total pulse is plus or minus 
the identity so that a could take any direction. This 
is reflected in the singuarities of '0 = 2TO7r, m G Z. 

At i = the current axis of rotation v and the 
global one a coincide. The former lies by construction 
in the xy-plane implying the inital conditions 



lim6'(t) 
t->o 



TT 

2' 



(34a) 



limthU) = 0, 

t-i-O 

lim0(t) 0(0). 



Note that the last two equations represent our arbi- 
trary choice which direction is the a;-direction for the 
spins and where the phase fl starts. In the numerical 
solutions below we will use f2(0) ~ 0. Inspecting the 
limit t — )■ one additionally finds 



dt 



t=o 



dt9 |t=o = 0. 



(35a) 
(35b) 



The derivative dti}} follows trivially from Eq. p3ap . 

We aim at solutions of Eqs. (|15ap and (jl5b[) both 
for TT and 7r/2 rotations and for pure dephasing fj [t) = 
(0,0,r] (t)). The corresponding vector n^(i) read 

Hfjit) = ry (t) Hy (t) 



-ay sin (jp) + (1 — cos (■0)) OxOz 



ax sin (V') 
cos (7/)) 



(1 - COS (0)) aya^ 
(1 - cos {%lj))al 



(36) 



For the sake of simplicity we have omitted the explicit 
time dependence of a and ^j. 




0.4 a 



Figure 5: (Color online) First order tt pulse suppressing pure 
dephasing by frequency modulation. Parameters given in Tab. 

For the phase ^{t) we use the Fourier series ansatz 



(34b) r!(t) = ^62„-isin 

n 

(34c) 



, 2-Knt , 



- &2n[ 



cos I 



(37) 
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with 



A*i,2 = Vz I '^yz (t) dt 



Figure 6: (Color online) First order 7r/2 pulse suppressing pure 
dephasing by frequency modulation. Parameters given in Tab. 

m 



1st order FM tt- and ^ pulses 

TT-pulse ^ -pulse 

~Vq 3.751157 Vo 4.928277 

bi 5i 

62 -1.090479 b2 -0.944852 

63 63 

64 -0.588913 64 -0.122088 

Table 1: Parameters for the first order frequency modulated 
(FM) pulses suppressing pure dephasing by satisfying Eqs. I I41II 
for compensating a classical bath. The coefficients can be com- 
pared to those derived in Ref. [30ll for a quantum bath. The 
dimensionless coefficients bn refer to the ansatz JSTjl. The am- 
plitudes Vo are given in units of 1 /rp . The odd coefficients are 
numerically within fO^'^'^. 



The value 9{Tp) is fixed by the fact, that the final axis 
of rotation has to be perpendicular to to rotate the 
spin by the full angle il){Tp). Thus one has 



(rp) = f (38) 



In first order, explicit computation yields 



TpHl = m^iax + ^ll,2^y + Ail,30-z (39) 



(40a) 
(40b) 
(40c) 



where the functions naj3{t) with q:,/3 G {x,?/, z} are 
given by the matrix elements of the rotation matrix 
Dai—'ip) in (|B.1[) in [Appendix B[ Thus, for first order 
pulses one has to achieve 



Aii,3 ^Vz '"'zz (Odt, 



(41) 



for i G {1,2,3}. Typical solutions for tt and tt/2 
pulses are shown in Figs. [5] and IH] The parameters 
are given in Tab. [TJ 

In second order we similarly obtain 



(42) 



with 



A*2,l 



A*2,2 



/^2,3 



Tp ti 





Tp t 





Tp ti 



{nyz,inzza - '^zz,i%z,2)dt2dti (43a) 
{nzz,inxz,2 - nxz,inzz,2)<^t2dti (43b) 
{nxz,inyz^2 - "-X2.i'^xz,2)dt2dti, (43c) 







using the shorthand Uap^k '■= n^p (tk) where k G 
{1,2}. Thus a second order FM pulse has to fulfill 
(|iT|) and 

= ^l2,^. (44) 

Solutions for second order pulses with frequency 
modulation are shown in Figs. [7] and |8] for tt pulses 
and in Figs. [9] and [TOl for tt/2 pulse. In Figs. [7] and |9] 
the phases ^l{t) arc shown while in Figs. 151 andlTUlthe 
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Figure 7: (Color online) Minimized second order FM tt pulse 
suppressing pure dephasing with Q{t) as in Eq. II37I I. The 
cocfBcients for this pulse are given in Tab. [2] 

corresponding amplitudes are plotted. These pulses 
are minimized in the following way. Since it is impor- 
tant to realize pulses with small amplitudes for prac- 
tical purposes, we studied sets of Fourier coefficients 
in (|37p with one coefficient more than necessary for 
the conditions (|4T|) and (|44| . We used this freedom 
to minimize the amplitude and carried this out for 
four different choices of the Fourier coefficients. The 
solution with the minimal amplitude out of the four 
minima is shown. The characteristics of the pulses 
are reported in Tabs. [2] and [3l respectively. We stress 
again that the classical case requires less coefficients 
than the quantum case. The maximal amplitudes of 
the classical pulses are smaller than the amplitudes 
for the quantum pulse, as can be seen comparing 
the values in the tables or the plots for and Vy. 
Moreover, the pulses suppressing classical noise have 
a simpler form. 

Obviously, the solutions are symmetric because all 
sine coefficients in Tab. [2] vanish to numerical accu- 
racy. This holds for the classical as well as for the 
quantum case. Deviations from this symmetric shape 
in the results of Ref. [s^ are due to the lower preci- 
sion of the numerics used in this preceding article. 

In addition to the tools employed already for 
AM pulses, finding frequency modulated pulses re- 
quires to solve a system of three coupled ordi- 



Figure 8: (Color online) Amplitudes for the same pulse as in 
Fig.El 

nary differential equations in each step of the mul- 
tidimensional root search. This was done by us- 
ing "gsLodciv2_system" with a stepping of the type 
"gsLodeiv2_step_rk4" . Thus a fourth order Runge- 
Kutta integration with adaptive step size governed 
by the double step method to keep the local absolute 
error estimate in the order of magnitude of 10~^^. 

5. Amplitude and Frequency Modulated 
Pulses 

Allowing for the modulation of amplitude and fre- 
quency leads to a humongous parameter space. Thus 
we restrict ourselves to illustrating that the frequency 
modulated pulses of the previous section can be mod- 
ified to account for smooth transients when the pulses 
are switched on and off. Such transients are generic 
in experimental realizations. 

To describe the transient region we define 

r /onW < i < T, 

fit) < I Ts<t<Tp-Ts 

[ 1 - .fon[t - (1 - r,)] Tp-Ts<t<Tp 

(45) 

where /on(i) := sin^ {nt/2Ts). The actual amplitude 
V{t) is given by Vof{t). Then Tg stands for the time 
needed to reach the amplitude Vq from zero ampli- 
tude or vice- versa, see Fig. [TTl 
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Minimized 2nd order tt FM pulses 



Minimized 2nd order 7r/2 FM pulses 









mi^iTitii iTi 








nnpi nfmn 


Vn 


a 1 9QnQ7 


•'0 


1 7071 1 5 


Vn 


7 405785 


Vn 


8 43541 4 


bi 





bi 





bi 


1.524556 


61 


-1.820216 


62 


-0.381075 


b2 


1.392956 


62 


-0.349899 


62 


-0.351972 


63 





b3 





63 


0.325909 


63 


0.030436 


64 


0.450018 


bi 


-0.705159 




0.411212 


64 


0.521648 


65 





65 





65 


0.690512 


65 


-0.555341 


b6 


-0.496673 


be 


0.133042 


66 


-0.510771 


be 


-0.387557 


67 





h 





67 


0.347745 


67 


0.451462 




-0.241963 


bs 


0.690594 


611 


0.019634 


68 


-0.193733 






b9 









69 


-0.161450 






bw 


-0.695501 






610 


-0.282067 






614 


0.472195 






614 


0.047116 



Tabic 2: Parameters of the minimized FM tt pulses suppressing 
pure dephasing by satisfying the first and second order Eqs. 
I I41II and II44I I for a classical or the corresponding conditions 
for a quantum dephasing bath in Rcf. [30l |. The dimensionless 
coefficients b„ refer to Eq. II37I I. The amplitudes Vb are given 
in units of l/rp . The odd coefficients are numerically within 
10-1". 



Table 3: Parameters of the minimized FM 7r/2 pulses sup- 
pressing pure dephasing by satisfying all first and second or- 
der Eqs. II41I I and II44I I for a classical bath or the corresponding 
conditions for a quantum bath in Ref . [30l | . The dimensionless 
coefficients 6,1 refer to Eq. II37I I. The amplitudes Vq are given 
in units of l/rp. 



In this section, we report the results for Tg = O.lTp. 
In [Appendix A[ results for other values of Tg arc in- 
cluded for comparison. The parametrization of Q{t) 
and a{t) are the same as in section [J] 

Examples for first order AM-fFM tt and it/2 pulses 
are plotted in Figs. [T2l and [T3l and the coefficients of 
the parametrization can be found in Tab. [H The 
maximum amplitude of the AM-I-FM pulses is lower 
than the amplitude of the continuous AM pulses and 
of the FM pulses derived for a quantum bath [2^. 
It is comparable to the amplitude of the well-known 

which is 77r/6 



SCORPSE TT pulses 
-1 



35, 3 



3.7 in 



units of T~ . In comparison to the pure FM pulses 
of first order shown in Figs. [5] and [6] with parameters 
given in Tab. [1] we clearly see that a bit larger ampli- 
tudes arc required for finite transients as was to be 
expected. 

Analogous considerations hold for second order 
pulses. They are plotted in Figs. [14] and [TSj Their 
coefficients are listed in Tab. \5[ 



6. Frequency Modulated Pulses for General 
Decoherence 

In section [4] on frequency modulated pulses we 
pointed out that one of their major assets is that they 
realize rotations about two independent spin axes. 
Hence FM pulses may correct not only dephasing 
without spin flips but also longitudinal relaxation in- 
cluding spin flips. We treat this situation here. Thus 
we deal with the Hamiltonian in ([5]) where all three 
components rja {a € {x, y, z}) of ff are present. 

We consider a bath with cylindrical symmetry 

First order tt and tt/2 AM-I-FM pulses 

TT Tt/2 



Vo 
61 
62 
63 
64 



4.232216 


-1.073059 


-0.233720 



Vo 
61 
62 
63 
64 



5.930552 


-0.506131 


0.053241 



Table 4: Parameters for the tt and n/2 AM-I-FM pulses sup- 
pressing pure dephasing with Tb = O.lTp. The coefficients bn 
refer to Eq. H37|l . The amplitudes Vb are given in units of 1/rp . 
The odd coefficients vanish numerically with in 10"^^. 
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Figure 9: (Color online) Minimized second order FM 7r/2 pulse Figure 10: (Color online) Amplitudes for the same pulse as in 
suppressing pure dcphasing with Q{t) as in Eq. I|37|l . The Fig.|9l 
coefficients for this pulse are given in Tab.|3l 

generic for an NMR experiment where the z-axis of 
the spin is distinguished by a large magnetic field and 

the system is rotationally invariant around this axis, cross correlations to conclude that the vanishing of 
Considering rotations by tt and ^ about cr^ one easily 
finds that the averages of rf fulfill 



Vx = Vy 







(46) 



Additionally, second order pulse require the autocor- 
relation 



TpH2 = (A*2,l+/^2,4)o"x + (A*2,2+M2.5)o-y + (M2.3+Ai2,6)a'z 

(48) 

requires the coefficients fi2,i to take the value zero 



(47) 



with a 6 {x, y, z}. We do not mention terms linear in 
\t\ because we aim at second order pulses at most. In 
principle, the cross correlations ga^p{t = 0) may also 
matter. But the cylindrical symmetry in combination 
with the antisymmetry of g^yi—t) ~ ^Qxyit) imply 
that all cross correlations ga^p{t = 0) vanish. In 
addition, = gxx(^) = SyyiO) is implied. 

On the basis of the expectation values ^ , the van- 
ishing of the first term, sec Eq. p9p . requires the van- 
ishing of the same integral equations as in the case of 
pure dephasing (|40p for frequency modulated pulses. 
Hence the same solutions result presented in the two 
previous sections. 

In second order, we use the variances and vanishing 



Tp ti 

fJ-2,1 = sl J J [nyx,inzx,2 + nyy^in:,y^2 



- {,nz^snyx,2 + n^y^inyy^2)]dt2dti (49a) 

■^P ti 

,2, = 4j l[n^^,n^.^, + n.y,n^y,2 



- ("-xx,i"-zx,2 + n^y sn2y^2)]dt2dti (49b) 

Tp ti 

1^2,3 = y J [nxx,inyx,2 + nxyS''T'yy^2 


- {nyx.inxx,2 + "-TO4"xy,2)]dt2dti (49c) 



13 



0.2 0.4 0.6 0.8 1 

t [T ] 

Figure 11; (Color online) A more realistic variation of the 
amplitude V{t) = Vof{t) according to Eq. I|45p for frequency 
modulated pulses. Note the transients for switching the pulse 
on and off. 



Tp ti 







4) 





(i2) 










nyz (i2)]dt2dti 

Tp ti 




(49d) 






4) 


j j[nzz{t.)n.z 




fe) 










n^, (t2)]dt2dti 

Tp ii 




(49e) 




{vl + 




j j [n-xz ih) 




(i2) 








ih] 


n-xz (i2)]dt2dti. 




(49f) 



Note that we have to require that each ^2,i vanishes 
even though they appear in pairs in front of the Pauh 
matrices in (|48p because we do not want to make 
assumptions of the relative size of and rf^ + si- 

The parametrization of Vl{t) and a{t) are chosen as 
before in sections |3] and [SJ Solutions to the condi- 
tions ^i^i = and /i2,i = are shown in Fig. 1161 for 
a IT pulse and in Fig. [T7] for a 7r/2 pulse. The pa- 
rameters of these pulses are reported in Tab. [51 The 
amplitudes of these pulses have been minimized as 
described before. 

Comparing the amplitudes in Tab. [5] with those for 



Figure 12: (Color online) AM+FM first order tt pulse sup- 
pressing pure dephasing for Ts = O.lrp. The phase f2(f) is 
parametrized as in Eq. I I37II . The coefficients of this pulse are 
given in Tab. 2] 

classical pulses in Tabs. [5] and |31 it is clear that the 
suppression of general decoherence requires higher 
amplitudes or longer pulses, respectively. But the 
increase in amplitude is not very large. For the tt 
pulse, the amplitude is increased by 12% and for the 
7r/2 pulse it is even lowered by 1%. This finding ap- 
pears contradictory at first sight because more condi- 
tions have to be fulfilled for general decoherence. The 
contradiction is resolved by the observation that we 
consider more coefficients in the construction of the 
pulse suppressing general decoherence than we do in 
the construction of the pulse suppressing pure de- 
phasing, see Tabs. [3] and [6l Of course, the 7r/2 pulse 
suppressing general decoherence fulfills also the con- 
ditions for the suppression of pure dephasing. It is 
remarkable that the additional conditions for general 
decoherence can be fulfilled at only moderate addi- 
tional effort. Hence one may realize a pulse which 
not only suppresses transversal decoherence but also 
longitudinal decoherence. 

Another intriguing possibility is to use the pulses 
found as replacement for XY4 and XY8 cycles 0, 0] ■ 
The main difi^erence between the net effect of XY4 to 
the TT pulse suppressing general decoherence in Fig. 
[TBI and Tab. |6] is that XY4 is designed to be a no- 
operation (NOOP) sequence while the tt pulse realizes 
a spin flip. But if the tt pulses is applied a second time 
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2nd order AM+FM tt and 7r/2 pulses 



0.3 




n 




7r/2 






9.076304 


Vo 


10.450781 


0.25 


h 





hi 





0.2 g 


b2 


-0.436689 


b2 


-0.123441 


b3 





hi 





0.15^ 


hi 


0.305937 




-0.130381 









b5 





0.1 


be 


-0.585209 
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-0.679511 
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Figure 13: (Color online) AM-I-FM first order it/2 pulse sup- 
pressing pure dcphasing for Tg = O.lrp. The phase Q{t) is 
parametrized as in Eq. l|37|l . The coefficients of this pulse are 
given in Tab. [4] 



in time-reversed order f2(T) — Q{Tp — t) it realizes 
again a second order tt pulses. Thus the effect of both 
pulses back-to-back is a 27r pulse which reduces to a 
mere phase factor —1 without effect on the actual 
spin state. Hence, such a composite pulse of length 

2Tn 



FM, forward 



FM. back ward 



(50) 



is symmetric and corresponds to an XY8 cycle in the 
sense that it suppresses spin dcphasing and spin flips; 
thus it suppresses general decoherence up to second 
order. The intriguing aspect is that this composite 
pulse is an always-on pulse and does not consist of 8 
individual pulses. Thus its amplitude is much lower 
than the amplitude needed in an XY8 cycle. For 
given pulse length one can reduce the total energy 
cx VqTp needed for the coherent control. For given 
maximal amplitude the cycle can be performed much 
more rapidly. This opens a promising route to more 
effective control calling for experimental verification. 

7. Conclusions 

Coherent control of quantum systems is a very ac- 
tive field of current research. The simplest quantum 
system to be controlled generally is a two-level sys- 
tem which can be seen as a spin S = 1/2. In par- 



Table 5: Parameters for the tt and it/2 AM-|-FM pulses sup- 
pressing pure dephasing by satisfying the first and second order 
Eqs. 1141 I I and II44I I for a classical bath. The coefficients b„ refer 
to Eq. II37I I. The amplitudes Vq are given in units of l/rp. The 
odd coefficients vanish numerically within 10^^. 



ticular, the coherent control of such spins is at the 
very heart of magnetic resonance techniques, where 
nuclear spins are manipulated, and of quantum infor- 
mation processing, where quantum bits are manipu- 
lated. Previously, pulses were designed to suppress 
the influence of a noisy quantum environment. Al- 
though this is the most general case it is cumbersome 
because many conditions have to be fulfilled and thus 
the necessary pulses are fairly complicated. They 
need relatively large amplitudes or they are relatively 
long. Both properties limit their practical relevance. 

In practice, however, very often the environmen- 
tal noise is dominated by classical fluctuations, for 
instance because it results from macroscopic degrees 
of freedom at relatively high temperature (room tem- 
perature) . This led us to study pulses which suppress 
classical noise. Such pulses arc subject to less condi- 
tions so that they can be simpler and of lower ampli- 
tudes or their duration Tp can be shorter. We studied 
such pulses analytically by a systematic expansion in 
Tp. Pulses with amplitude modulation, which rotate 
the spin about a single fixed axis, and pulses with 
frequency modulation, which rotate the spin about a 
continuously varying axis in the xy-plane of the spin 
directions, are considered. 

A purely dephasing bath is relevant if the energy 
splitting between the two levels is so large that the 
rotating frame approximation is applicable. Exper- 
imentally, the dephasing dominates if T2 is much 
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Figure 14: (Color online) Minimized second order AM+FM 
TT pulse suppressing pure dephasing with Q{t) as in Eq. I I37II . 
The coefficients are given in Tab. [3] 



Figure 15: (Color online) Minimized second order AM+FM 
tt/2 pulse suppressing pure dephasing with Q{t) as in Eq. II37I I. 
The coefficients are given in Tab. [S] 



smaller than Ti . For this situation we presented first 
and second order pulses based on amplitude or on fre- 
quency modulation, respectively. We explicitly pre- 
sented pulse shapes for tt and tt/2 pulses which either 
flip the spin between up and down or which rotate 
them by 90° from the z-axis. In all pulses suppressing 
classical noise we confirmed that their amplitudes are 
lower than for their previously known quantum coun- 
terparts. The only exception are first order pulses. In 
this order the quantum character is not manifest in 
the conditions so that they are identical for classical 
and for quantum baths. 

Furthermore, we showed that combinations of 
amplitude and frequency modulation can also be 
treated. In particular, finite transients in the switch- 
ing processes of amplitudes can be accounted for. 
Thus imperfections in the switching can be consid- 
ered and they do not pose a conceptual problem. 

Intriguingly, we could furthermore establish the ex- 
istence of pulses which suppress classical general de- 
coherence. This means that not only transversal de- 
phasing but also longitudinal spin relaxation relying 
on spin flips can be suppressed. Amplitude modula- 
tion is not sufficient to this end because one needs at 
least two independent axes of rotation to suppress all 
kinds of spin errors. But frequency modulated pulses 
can do the job. 

We found that the necessary amplitudes are at 



worst only moderately larger than for the suppres- 
sion of pure dephasing alone. Thereby, we propose 
a single shaped pulse which has similar properties 
as cycles of pulses. In particular, a composite pulse 
of two TT pulses, which suppress general decoherence, 
can replace the well-known XY8 cycle. The asset of 
the composite pulse we are advocating here is that it 
is an always-on pulse. Hence the required amplitudes 
are much lower than those required for a cycle of very 
short pulses. 

Of course, further research is called for. On the 
experimental side, studies of the performance of the 
proposed pulses are called for. A key question is 
whether the proposed shapes can be realized reli- 
ably enough to reach the predicted suppression of the 
noisy baths. 

Theoretically, the question of the robustness of the 
proposed pulses towards imperfections in their real- 
izations deserves investigations. For instance, simu- 
lations of the pulses in various baths are called for 
to guide experiment. In parallel, we believe that the 
necessary amplitudes can still be reduced by further 
minimization. 

Finally, we point out that analogous expansions 
for coupled two-level systems have hardly attracted 
attention so far in spite of their relevance of two-qubit 
gates in quantum information processing. 



16 




Figure 16: (Color online) Minimized second order FM tt pulse 
suppressing general decoherence with Q{t) as in Eq. I I37II . The 
coefficients are given in Tab. |6] 



Figure 17: (Color online) Minimized second order FM n/2 
pulse suppressing general decoherence with Q{t) as in Eq. I|37|l . 
The coefficients are given in Tab. |6] 
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Appendix A. Pulses with Amplitude and Fre- 
quency Modulation 

In section \5\ we showed an example of two FM 
pulses with angle tt and 7r/2 and finite transients of 
amplitude for switching on and off. Each transient 
takes a fraction Tg/rp of the total time Tp of the pulse. 
Here we provide further tt pulses with shorter tran- 
sient time in Tab. IA.7I in first order and in Tab. IA.8I 
in second order. The parameters for the first order 
pulses are to be compared to those reported in Tab. ID 
Remarkably, the amplitude does not increase mono- 
tonically on increasing Tg/rp in Tab. IA.7I 

The second order parameters illustrate that ampli- 
tude and frequency modulation can also be combined 
in second order. The parameters can be compared to 
those of pure FM pulses given in Tab. [2] 

Appendix B. Rotation Matrix 

For the derivation of the matrix Z?a(— -0) we re- 
fer the reader to Ref. [l^. We obtain the matrix 
()B.ip below, where the time dependencies of ip{t) 
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2nd order FM n and n/2 pulses for general decoherence 
TT 7r/2 



Vo 


9 


.079728 




7, 


,297361 


h 


1 


.818085 


61 




1.793195 


b2 





.514273 




0, 


.223583 


b3 




0.231238 




0, 


.221590 


bi 




0.220323 


bi 


0, 


.324311 


65 


0, 


.014857 






0.579783 


be 


0, 


.508720 




0, 


.272144 


br 




0.439837 


67 


0, 


.507358 


bs 




0.816150 






0.119786 


b9 




0.332547 


69 




0.011429 


bio 




0.846412 


610 


0, 


.069581 


bii 




0.249481 


613 


0, 


.219071 



Table 6: Parameters for the minimized FM vr and 7r/2 pulses 
suppressing general decoherence. The coefficients b„ refer to 
Eq. 13711 . The amplitudes Vq are given in units of l/vp. 



1st order AM+FM Trpulses 

Ts = O.OlTp = O.OOlTp 

Vo 3.907279 Vo 4.016102 

61 5i 

62 -0.892324 62 -0.791296 

63 &3 

64 -0.196455 64 -0.040110 

Table A. 7: Parameters of AM+FM first order tt suppressing 
pure dephasing for two switching times Ts. The dimensionless 
coefficients bn refer to Eq. I)37|l. The amplitudes Vb are given 
in units of l/rp. The odd coefficients vanish within fO~^^. 



2nd order AM+FM Trpulses 



0.1r„ 



O.OIto 



Vo 9.076304 

62 -0.436689 

64 0.305937 

bo -0.585209 



Vo 
62 
64 
bo 



8.486171 
-0.309163 
0.507966 
-0.437161 



Table A. 8: Parameters of AM+FM second order vr suppressing 
pure dephasing for two switching times Tg. The dimensionless 
coefficients bn refer to Eq. JST}. The amplitudes Vb are given 
in units of l/rp. The odd coefficients vanish within fO"''. 
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008-0 + (1 — cosip)a'^ 
Dai-i') = { -a^ sin-0 + (1 - COS V')a£!;ay 
ay sin ■(/) + (1 — cos ^)axaz 



az sin ip - 



(1 
(1 



cosip)a'i 



-ax sin Ip + (l — cos 'ijj)ayaz 



—tty sin '0 + (1 — cos ip)axaz 
ax sin + (1 — cos '4>)ayaz 
cos^p + (1 — cosijj)al 



(B.l) 
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